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We measure the rate of D'^-D'' mixing with the observable ycp = (tktt/tkk) — 1, where tkk and tktt 
are respectively the mean lifetimes of CP-even D° — > K'^K~ and CP-mixed D'^ — > K^tv^ decays, using 
a data sample of 384 fb~^ collected by the BaBAR detector at the SLAC PEP-II asymmetric-energy B Factory. 
From a sample of and D" decays where the inital flavor of the decaying meson is not determined, we obtain 
ycp = [1.12 ± 0.26(stat) ± 0.22(syst)]%, which excludes the no-mixing hypothesis at 3.3a, including both 
statistical and systematic uncertainties. This result is in good agreement with a previous BaBAR measurement 
of ycp obtained from a sample of D*"*" D^tt"*" events, where the decays to K~-jt^, K~ , and 
TT+TT^, which is disjoint with the untagged D'^ events used here. Combining the two results taking into account 
statistical and systematic uncertainties, where the systematic uncertainties are assumed to be 100% correlated, 
we find ycp = [1.16 ± 0.22(stat) ± 0.18(syst)]%, which excludes the no-mixing hypothesis at 4.1(t. 
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Several recent results Uj, |2|, |3l |4l] show evidence for mixing 
in the D^-D^ system consistent with predictions of possible 



These results also 
[ni|ll[ll[ll,and 



Standard Model contributions iSSItIJE. 
constrain many new physics models jlOl | 
increasingly precise D^-D'^ mixing measurements will pro- 
vide even stronger constraints. One manifestation of D^-D^ 
mixing is differing decay time distributions for decays to 
different CP eigenstates [15]. We present here a measurement 
of this lifetime difference using a sample of and decays 
in which the initial flavor of the decaying meson is unknown. 

Assuming CP conservation in mixing, the two neutral D 
mass eigenstates \Di) and \D2) can be represented as 



\D2) = 



p\D^) 



q\D^) 



(1) 



where \p\^ + \q\^ = 1. The rate of D'^-D'^ mixing can be 
characterized by the parameters x = Am/F and y = Ar/2r, 
where Am = nii — 7712 and AF = Fi — F2 are respectively 
the differences between the mass and width eigenvalues of the 
states in Eq. ([T]i, and F = (Fi + F2)/2 is the average width. 
If either x or y is non-zero, mixing will occur, altering the 
decay time distribution of D° and D'^ mesons decaying into 
final states of specific CP lfl6ll . 

In the limit of small mixing, and no CP violation in mixing 
or in the interference between mixing and decay (assumptions 
which are consistent with current experimental results), the 
mean lifetimes of decays to a CP eigenstate of a sample of D*^ 
{r^f^ ) and (r^ ), along with the mean lifetime of decays 
to a state of indefinite CP {tkt,), can be combined into the 
quantity 



ycp 



1, 



(2) 



where (r,,, 



D" 
hh 



/^")/2. Noting that the untagged 



A'^7r+ II17II final state is a mixture of Cabbibo-favored and 
doubly Cabbibo-suppressed Z?" and decays with a purely 
exponential lifetime distribution, along with a very small ad- 
mixture of mixed D° decays, an analogous expression also 
holds for {tktj). Given the current experimental evidence in- 
dicating a small mixing rate, the lifetime distribution for all 
hh and Ktt final states is exponential to a good approxima- 
tion. If ycp is zero there is no D'^-D^ mixing attributable to 
a width difference, although mixing caused by a mass differ- 
ence may be present. In the limit of no direct CP violation, 

ycp = y- 

We measure the mean lifetime in the decay modes 
K~TT^ and K'^K~ , where the initial flavor of the decay- 
ing D° is not identified (the untagged sample). This sam- 
ple excludes mesons which can be reconstructed as part 
of D*^ Z?°7r+ decays, as these decays (the tagged sam- 
ple) are the subject of an earlier BaB/^ analysis ( ISQ whose 
results are combined with those of the current analysis. To 
avoid potential bias, we finalized our data selection criteria, 
fitting methodology, sources of possible systematic uncertain- 
ties to be examined, and method of calculating statistical lim- 
its for the current untagged analysis alone and in combination 



with the tagged analysis, prior to examining the mixing results 
from the untagged data. In general, systematic uncertainties 
related to the reconstruction of signal events cancel in the life- 
time ratio. However, uncertainties related to the somewhat 
differing backgrounds present in the /\ ^7r+ and /\ +/\ ^ final 
states lead to larger systematic uncertainties in the untagged 
analysis compared to those of the tagged analysis, which has 
much higher signal purity . 

We use 384 fb~^ of e+e^ colliding-beam data recorded 
at, and slightly below, the T(4S') resonance (center-of-mass 
[CM] energy ^ ^ 10.6 GcV) with the BaR4R detec- 
tor [[Q at the SLAC National Accelerator Laboratory PEP-II 
asymmetric-energy B Factory. Candidate signal decays 
are reconstructed in the final states A'^7r+ and K^K^ . The 
selection of events and reconstruction of D° signal candi- 
dates closely follows that of our previous tagged analysis ifisl] . 
We require A' + and tt^ candidates to satisfy particle identi- 
fication criteria based on cLE / da; ionization energy loss and 
Cherenkov angle measurements. We fit oppositely charged 
pairs of these candidates with appropriate mass hypotheses to 
a common vertex to form a candidate. The decay time 
t of each candidate with invariant mass within the range 
1.80 — 1.93 GeV/c^, along with its estimated uncertainty at, 
is determined from a combined fit to the D'^ production and 
decay vertices, with a constraint that the production point be 
consistent with the e+e^ interaction region as determined on 
an event-by-event basis. We retain only candidates with a 
X^-based probabiHty for the fit P{x^) > 0.1%, and with 
-2 < < < 4 ps and at < 0.5 ps. 

We further require the helicity angle 9h, defined as the 
angle between the positively charged track in the rest 
frame and the D'^ direction in the laboratory frame, to sat- 
isfy |cos6'/f| < 0.7, which aids in the rejection of purely 
combinatorial background events. Contributions from true D'^ 
mesons produced in B meson decay are reduced to a negligi- 
ble amount by rejecting D'^ candidates with momentum in the 
e+e^ CM frame less than 2.5 GcV/c. For events with multi- 
ple candidates sharing one or more tracks, we retain only the 
candidate with the highest P(x^). The fraction of events with 
multiple signal candidates is ~ 0.05% for the K^K^ final 
state, and ~ 0.3% for K^n^ . 

The invariant mass distributions for the final D'^ A'^7r+ 
and D'^ K^K^ samples are shown in Fig.[I] For the life- 
time fits, we use only events within ±10MeV/c^ of the 
signal peak 1.8545 < M^o < 1.8745 GeV/c^ (the lifetime 
fit mass region). The K^tt^ and AT+AT^ signal yields within 
this region and their purity are given in Table|T] Events within 
the mass sideband regions 1.81 < Mj^o < 1.83GeV/c^ 
and 1.90 < Md« < 1.92GeV/c^ are used to determine the 
combinatorial background decay time distribution within the 
lifetime fit mass region. In addition to purely combinatorial 
backgrounds, there are small background contributions from 
decays of non-signal charm parents where two of the decay 
products are selected as the daughters of a signal decay and 
subsequently pass the final event selection. These misrecon- 
structed charm backgrounds are accounted for using simu- 
lated events. Their contribution is ^ 0.7% (^ 3.8%) of the 
total number of background events in the K~tt~^ {K^K~) 
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FIG. 1: (a) D° K--K+ and (b) D° -> /-C+if - invariant mass dis- 
tribution with the data (points), total fit (line) and background contri- 
bution (solid) overlaid. The innermost dashed lines on either side of 
the signal peaks delimit the lifetime fit mass region, with lower and 
upper mass sidebands shown on either side. 



TABLE I: D" K-tv+ and D° K+ R- signal yield and purity 
in the lifetime fit mass region. 



Sample Signal Yield (x 10^) Purity (%) 



K--K+ 
K+K- 



2710.2 ±3.4 
263.6 ± 1.0 



94.2 
80.9 



signal region. 

The mean Z?" lifetime is determined from a fit essentially 
identical to the one performed in the previous tagged analy- 
sis ifisll . using the reconstructed decay time t and the decay 
time uncertainty at for events within the lifetime fit mass re- 
gion. Three categories of events are accounted for in the life- 
time fit: signal decays, combinatorial background, and misre- 
constructed charm events. 

The decay time distribution of signal events is described 
by an exponential convolved with a resolution function which 
is taken as the sum of three Gaussian functions with widths 
proportional to at- The functional form of this probability 
density function (PDF) for signal events is 

TZx{t,cTt;Tx) = ft3'D{t,(Tt; Sxs3,to,Tx) 

+ - fta) ft2'D{t,at;SxS2,to,Tx) (3) 

+ {I - ft2)T^it,o-t;Sxsi,to,Tx) , 

where fu (with i = 1 ... 3) parameterizes the contribution 
of each individual resolution function, Si is a scaling factor 
associated with each Gaussian, tx (where X = Kn, KK) is 
the lifetime parameter determined by the fit, is an offset to 
the mean of the resolution function, and where 



'D(t,at;s,to,T) = 

/ exp(-itruc/T)exp 



(t-ttr„c+to)" 



(4) 



with normalization coefficient Ccr^ . Up to an overall scale fac- 
tor in the width, the resolution function is identical for both 



final states. We account for a small (^ 1%) difference in the 
K~TT'^ and K^K^ resolution function width using an addi- 
tional fixed scale factor Sx- The value of Skk is determined 
from the data, with Sk-k fixed to 1 .0. Possible biases result- 
ing from this assumption are included as part of the study of 
systematic uncertainties. All other resolution function param- 
eters are shared among the two modes, and all parameters are 
allowed to vary in a simultaneous extended unbinned maxi- 
mum likelihood fit to both final states. 

The decay time distribution of the combinatorial back- 
ground is described by a sum of two Gaussians and a modified 
Gaussian with a power-law tail to account for a small num- 
ber of events with large reconstructed lifetimes. The widths 
of these Gaussians are not scaled using event-by-event uncer- 
tainties. Events in the lower and upper A'^7r+ (A'+ A'^) mass 
sidebands are fit separately, and a weighted average of the re- 
sults of these fits is used to parameterize the PDF for K^n^ 
(K^K^) combinatorial events in the lifetime fit mass region. 

Misreconstructed charm background events have one or 
more of the charm decay products either not reconstructed 
or reconstructed with the wrong particle hypothesis. In the 
K-TT+ (K+K-) final state, - 60% (- 95%) of these events 
are from true decays, with the balance coming from 
charged D and charm baryon decays. The charm background 
is long-lived and is described using an exponential convolved 
with a resolution function consisting of two Gaussians with a 
shared mean and widths that depend on at- Because the num- 
ber of these events in the K~'k'^ {K~^K~) sample is small 
relative to the total background, an effective lifetime distribu- 
tion taken from simulated events and summed over all A'^7r+ 
(AT+A'^) charm backgrounds is used in the A'^7r+ (A'+ AT^) 
lifetime fit. 

Since the lifetime fit PDFs depend on the event-by-event 
decay time uncertainty, PDFs describing the distribution of 
decay time uncertainties for each of the event classes are re- 
quired to avoid bias in the likelihood estimator used in the 
data fit ll20ll . We extract these distributions directly from the 
data. For combinatorial events, the distribution of decay time 
uncertainties is taken from a weighted average of the distri- 
butions extracted from the lower and upper mass sidebands. 
The decay time uncertainty distribution for signal events is 
obtained by subtracting the combinatorial background uncer- 
tainty distribution from the uncertainty distribution of all (i.e., 
background plus signal) candidates present in the lifetime fit 
mass region. The signal distribution is also used for the rel- 
atively small number of misreconstructed charm background 
events. 

The results of the lifetime fits are shown in Figs. |2] and [3] 
along with a plot of the point-by-point residuals for each 
fit normalized by the statistical uncertainty associated with 
a data point. We find the 1?° K^tt^ mean lifetime 
TK-K = 410.39 ± 0.38(stat) fs and the D° ^ A'+A'" mean 
lifetime tkk = 405.85 ± l.OO(stat) fs, yielding ycp = 
[1.12 ± 0.26(stat)]%. The statistical significance of this mix- 
ing result without taking into account systematic uncertainties 
is 4.3(7. This untagged result is in good agreement with our 
previous tagged analysis ifisll . When the two results are com- 
bined, we find ycp = [1.16 ± 0.22(stat)]%, a result with a 
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FIG. 2: D" K~Tv'^ decay time distribution with tlie data (points), 
total lifetime fit (line), combinatorial background (gray) and charm 
background (black) contributions overlaid. 
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FIG. 3: -D° —> K~ decay time distribution with the data 
(points), total lifetime fit (line), combinatorial background (gray) and 
charm background (black) contributions overlaid. 



statistical significance of 5.3cr, excluding any systematic un- 
certainties. 

Numerous cross-checks have been performed to assure the 
unbiased nature of the fit model and to validate the assump- 
tions used in its construction. We have performed fits to 
datasets composed of fully simulated signal and background 
events in the proportions seen in the actual data, and find no 
bias in the measurement of individual tk-k and tkk lifetimes 
for simulated signal events generated at 411.6 fs (very near the 
nominal Z?" K~tt~^ lifetime value lfT6ll ). or for a lifetime 
value ^ 10% greater than this for _D" K+K~ . We ad- 
ditionally find no significant variations in the reconstruction 
efficiency for signal decays as a function of the true decay 
time. 

Many of the systematic uncertainties associated with the in- 
dividual lifetime measurements cancel to a great extent in the 
ratio of lifetimes. We consider as possible sources of system- 
atic uncertainty: variations in the signal and background fit 
models, changes to the event selection, and detector effects 
that might introduce biases in the lifetime measurements. 

We test the assumption of a shared signal resolution model 
by separately fitting each mode using completely independent 
resolution functions, and assign as a systematic uncertainty 
the magnitude of the change | Aj/cp | in ycp relative to the re- 
sult of the nominal fit. We additionally perform the nominal 
fit using a double Gaussian signal resolution model, and simi- 
larly assign a systematic uncertainty. The total uncertainty as- 
sociated with the choice of signal resolution model is 0.016%. 

To estimate possible biases correlated with the extent and 



position of the lifetime fit mass region, the size of the mass 
window is varied by ±2 and ±5 McV/ without changing the 
mass region center, and the center is shifted by ±0.5 McV/c^ 
while retaining the nominal 20 MeV/c^ width. The total sys- 
tematic uncertainty obtained from variations in the lifetime fit 
mass window is 0.110%. 

The modeling of the misreconstructed charm background 
is taken from simulated events, and we vary the expected con- 
tribution from these events by ±15%(±5%) for the K~'k'^ 
{K^K^) final state. These bounds are conservatively as- 
signed based on the results of other BaRar charm analyses in 
which the background modes here are fully reconstructed, and 
in which data and simulated event yields are found to agree 
within a few percent. We additionally vary the effective life- 
time used in the charm background lifetime fit PDFs by the 
same percentages, which corresponds to > ±2(t in the statis- 
tical uncertainty given the number of simulated events used. 
The largest jAj/cpl value within each of these two classes of 
variations is assigned as a systematic uncertainty, 0.0585% 
for the normalization variations and 0.0624% for the effective 
lifetime variations, which are then added in quadrature. 

We account for a possible bias associated with obtaining 
the combinatorial lifetime PDF in the lifetime fit mass region 
from data in the lower and upper mass sidebands by fluctu- 
ating the PDF parameters taking into account the correlations 
and statistical errors resulting from the sideband fits. We con- 
struct 100 PDF variations for each of the lower and upper side- 
bands for each of the final states, and then perform the nom- 
inal lifetime fit using each variation. We separately compute 
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TABLE II: Systematic uncertainties. 



Uncertainty Source 


\Aycp\ (%) 


Signal resolution model 


0.016 


Mass window 


0.110 


Misreconstructed Charm model 


O.OoD 


Combinatorial PDF 


0.115 


at selection 


0.069 


Overlap candidate selection 


0.017 


Detector effects 


0.093 


Total 


0.216 



the RMS of the 100 Aycp values associated with each of the 
four sets of variations, and assign the largest RMS, 0.115%, 
as a systematic uncertainty. 

We evaluate systematic uncertainties associated with the se- 
lection of the final dataset by individually varying the selec- 
tion criteria. We change the maximum allowed decay time 
uncertainty by ±0.1 ps, and assign the largest |Aj/cp| value, 
0.069%, as a systematic uncertainty. We vary the way in 
which signal candidates that share tracks with other signal 
candidates are selected by removing all overlapping candi- 
dates, and separately by also retaining all such candidates, 
and again take the larger of the resulting two jAj/cpl values, 
0.017%, as a systematic uncertainty. 

We account for possible detector effects which might bias 
the lifetime ratio by using several different detector configura- 
tions to re-reconstruct simulated event samples with statistics 
greater than the actual data for each configuration. These con- 
figurations include vertex detector misalignments, along with 
boost and beamspot variations, whose extent is based on resid- 
ual uncertainties in studies of mu-pair and cosmic events. The 
misalignment configurations introduce changes of up to 4 fs 
in both KK and Ktt lifetimes, as well as changes in the off- 
set parameter to of up to 5 fs. Since the same simulated event 
sample is reconstructed for each set of detector configuration, 
the variations are dominated by systematic effects. The total 
systematic uncertainty arising from this source is 0.093%. 

Table shows the contribution from each source of sys- 
tematic uncertainty given above. The total is calculated as the 
sum in quadrature of each of the individual items. In addition 
to the contributions quantified in the table, we also look for 
possible biases by fitting the data separated in: several differ- 



ent data-taking periods; several different azimuthal and polar 
angle bins in the laboratory frame for the D'^ candidate; sev- 
eral bins of the opening angle in the laboratory frame between 
the two D° daughters; several bins of the D'^ helicity angle; 
and several bins of the D'-^ momentum in the CM frame. We 
observed no significant biases in any of these cases. 

In our previously published tagged analysis ifisll . we com- 
bined the tagged result with the result of an untagged BaBar 
analysis done using a much smaller dataset ll2lll . and this pre- 
vious untagged result is superseded by the result here, which 
is ycp (untagged) = [1.12±0.26(stat)±0.22(syst)]%, which 
excludes the no-mixing hypothesis at 3.3ct, including both sta- 
tistical and systematic uncertainties. Our previous tagged re- 
suhiHis ycp(tagged) = [1.24±0.39(stat)±0.13(syst)]%. 
These results contain no events in common, and are thus sta- 
tistically uncorrected by construction. However, the degree of 
correlation in the systematic uncertainties is substantial, and 
we conservatively assume a 100% correlation in the system- 
atics shared between the two analyses. Combining the tagged 
and untagged results taking into account both statistical and 
systematic uncertainties 1221, we find ycp (correlated) = 
[1.16 ± 0.22(stat) ± 0.18(syst)]%. Summing statistical and 
systematic uncertainties in quadrature, the significance of this 
measurement is 4.1cr. 
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